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ABSTRACT 

The  stretching  and  break-up  of  a  viscoelastic  filament  pulled  by  a 
constant  weight  is  studied  numerically  by  a  finite  difference  method.  Hie 
results  show  the  following  tendencies: 

1.  Newtonian  filaments,  even  in  the  absence  of  surface  tension,  show  a  rapid 

M.1%*-  y  oX 

increase  in  elongation  at  one  particular  point  (they  "break"  there ) . 

2.  The  addition  of  a  viscoelastic  polymer  prevents  or  at  least  delays  the 
break-up,  even  if  it  makes  only  a  small  difference  to  shear  viscosity. 

3.  Surface  tension  accelerates  break-up,  but  even  in  the  presence  of  surface 
tension  elasticity  has  a  stabilizing  effect. 
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SIGNIFICANCE  AND  EXPLANATION 


Elongations 1  flows  of  polymeric  fluids  have  received  much  attention  both 
because  of  their  usefulness  in  distinguishing  reasonable  constitutive  laws 
from  bad  ones  and  because  of  their  practical  importance.  A  question  which  if 
of  great  practical  significance,  but  is  poorly  understood,  is  if  and  when  thin 
filaments  or  sheets  of  a  liquid  will  break-up.  This  problem  arises  e.g.  in 
the  formation  of  fibers,  in  lubrication  and  in  problems  involving  atomization 
into  droplets.  Empirically,  it  is  well  known  that  viscoelastic  fluids  are 
much  more  stable  against  break-up  than  Newtonian  fluids  with  comparable 
viscosity,  but  an  adequate  theoretical  explanation  has  so  far  not  been  given. 


In  this  paper,  we  examine  a  viscoelastic  filament,  which  is  attached  at 
the  upper  end  and  pulled  by  a  constant  weight  on  the  lower  end.  This  problem 
is  motivated  by  an  experiment  of  J.  Matta.  The  polymer  is  assumed  to  satisfy 
a  particular  constitutive  law  known  as  the  "rubberlike  liquid".  The  equations 
for  this  problem  are  formulated  and  solved  by  a  finite  difference  method.  The 
results  show  that  Newtonian  filaments  break  up  much  more  easily  than  visco¬ 
elastic  filaments.  This  is  the  case  whether  or  not  surface  tension  is 
present.  As  expected,  surface  tension  accelerates  the  break-up,  but  even  with 


large  surface  tension,  the  viscoelastic  filament  is  still  more  stable  than  the 
Newtonian  one. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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A  FINITE  DIFFERENCE  STUDY  OF  THE  STRETCHING  AND 
BREAK-UP  OF  FILAMENTS  OF  POLYMER  SOLUTIONS 


Peter  A.  Markovich*  and  Michael  Renardy** 

1 .  Introduction 

It  is  well  known  empirically  [f],  [2J ,  [7],  [9],  [10],  [13]  that 
filaments  of  viscoelastic  fluids  are  much  more  stable  against  break-up  than 
filaments  of  Newtonian  fluids  with  comparable  viscosity.  However,  so  far  this 
has  not  been  well  understood  theoretically.  Linear  stability  analysis  of 
uniform  jets  shows  that  in  fact  capillary  instabilities  should  grow  faster  on 
a  viscoelastic  jet  than  on  a  Newtonian  jet  [1J,  [2],  {73,  [9].  This  has  led 
to  the  conclusion  that  the  breakup  of  jets  or  filaments  must  be  governed  by 
nonlinear  effects.  However,  we  are  not  aware  of  any  precise  calculations  in 
the  nonlinear  regime. 

In  this  paper  we  study  a  problem  that  is  motivated  by  an  experiment 
designed  by  J.  Matta  [8]  for  the  measurement  of  elongational  flow  properties 
in  dilute  polymer  solutions.  In  this  experiment,  the  fluid  is  slowly  extruded 
from  a  vertical  nozzle  and  forms  a  drop  at  the  nozzle  tip.  As  this  drop 
reaches  a  critical  size,  it  begins  to  fall  and  drags  a  filament  behind  itself. 

He  consider  the  following  idealized  model  problem:  A  filament  is  at  rest 
up  to  time  t  *  0.  For  t  >  0,  one  end  of  the  filament  Is  attached  (at  the 
"nozzle"),  the  other  is  pulled  by  a  constant  weight  (the  "drop").  The  visco- 
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elastic  fluid  in  the  filament  is  assumed  to  obey  the  constitutive  law  for  a 
-rubberlike  liquid-  [3],  [4],  [12]. 

We  use  the  following  notation:  x  denotes  the  position  of  fluid 
particles  in  a  reference  configuration,  which  may  or  may  not  agree  with  the 
initial  configuration.  It  is  assumed  that,  in  the  reference  configuration, 
the  filament  has  uniform  radius  6  and  length  2L,  i.e.  we  have 
x  e  .  By  u(x,t)  we  denote  the  spatial  position  of  the  particle  x  at 

time  t.  p  denotes  the  density  of  the  fluid,  n  the  Newtonian  contribution 
to  the  viscosity,  a  the  memory  function  of  the  rubberlike  liquid  law,  a 
the  surface  tension  coefficient,  K  the  mass  of  the  weight,  and  g  the 
gravitational  acceleration.  The  problem  is  then  described  by  the  following 


(1.3)  u(x,t)  -  u0(x)  ,  t  <  0 

and  boundary  conditions 


(1.4)  u(-L,t)  *  «q <  —I. )  ,  t  >  0 

(1.5)  Mutt(L,t)  -  Mg  -  v«2s[ux] 

X*L 

These  equations  are  discretized  by  a  finite  difference  method,  which  is  an 


adaptation  of  a  scheme  used  earlier  on  a  similar  problem  [6] .  This  method  is 


described  in  section  2 


Zn  section  3  we  present  numerical  results.  Particular  emphasis  is  given 
to  the  question  oft  if  and  when  \i*  shows  catastrophic  growth.  We  shall  see 
that,  if  the  initial  thickness  is  uniform  (uq(x)  ■  x),  this  always  happens 
at  the  attached  end.  The  results  show  that  elasticity  retards  the  growth  of 
ux,  while  surface  tension  accelerates  it.  We  give  a  heuristic  explanation  of 
what  we  believe  to  be  the  mechanism  for  these  effects. 


2.  The  discretisation  scheme 

For  the  discretization  of  (1.1),  we  use  an  implicit  Euler  scheme  in  time 
and  sysnMtric  difference  approximations  in  space.  This  yields  a  stable  scheme 
which  is  first  order  accurate  in  time  and  second  order  accurate  in  space 
[6].  To  describe  the  scheme,  we  need  to  introduce  some  notations. 

The  grid  points  are  tn  ■  nk,  x^  **  ih  (i  «  -N,-N+1 , . . . ,N) ,  and 
approximations  to  u(xi#tn)  are  denoted  by  u”.  We  use  the  following 
difference  quotients: 


(2.1) 
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The  discretization  of  (1.1)  is  then  given  by 
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(2.2) 
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Of  course,  A+  and  A-  are  interchangeable.  For  j  <  0,  u^  is  simply 

Uq ( x j ) .  if  a  is  a  finite  sume  of  exponentials,  as  was  always  assumed  in  our 

calculations,  then  the  sums  can  be  evaluated  recursively. 

Equation  (2.2)  has  to  be  supplemented  by  boundary  conditions.  These  are 
n+1 


u 


-N 


u0(-D 


(2.3) 
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Here  A  is  a  second  order  accurate  approximation  to 

n 

*  n  V 

(2.4) 
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At  each  time  step,  (2.2),  (2.3)  yields  a  nonlinear  system  of  equations  for  the 


u”+1,  i  "  -N,...,N.  This  system  is  solved  by  Newton's  method. 


3.  Numerical  results 

In  order  to  test  our  program,  we  choose  the  following  parameter  values: 

L  «  1,  P  ”  1,  h  ■  10,  a(t)  “  e  t ,  y  “  1,  — ~  *  1,  g  *  1.  Me  used  the  initial 

3  *6* 

x 

condition  u.(x)  •  x  +  —  ,  and  we  added  inhomogeneous  terras  to  (1.1),  (1.4) 

x3 

and  (1.5)  in  such  a  way  that  u(x)  “  x  +  —  cosh  t  becomes  an  exact  solution. 

6 

Me  obtained  the  following  results  at  t  *  5: 
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h-0.05 

k-0.005 

a 

"x 

X  ■  -1 

-13.27 

38.08 

-0.8 

-7.H 

24.89 

-0.6 

-3.23 

14.43 

-0.4 

-1.15 

6.95 

-0.2 

-0.27 

2.48 

0 

0.03 

1.05 

0.2 

0.33 

2.48 

0.4 

1.21 

6.95 

0.6 

3.30 

14.47 

0.8 

7.19 

25.04 

1 

13.51 

38.67 

«xx 

-72.78 

-59.58 

-44.91 

-29.89 

-14.70 

0.00 

14.71 

29.97 

45.20 

60.58 

77.00 

h”0.0S 

k-0.001 

u 

«x 

-13.37 

38.13 

Eh 

m 

-1.18 

6.93 

-0.30 

2.48 

0.00 

1.04 

0.30 

2.48 

1.18 

6.93 

3.26 

14.38 

7.12 

24.84 

13.38 

38.19 

uxx 

-74.30 

B 

-44.77 

-29.78 

-14.66 

0.00 

14.66 

29.78 

44.79 

59.71 

74.76 

exact 

u 

ux 

-13.37 

38.10 

-7.13 

24.75 

-3.27 

14.36 

-1.19 

6.94 

-0.30 

2.48 

0 

1.0 

0.30 

2.48 

1.19 

6.94 

3.27 

14.36 

7.13 
24. 7S 

13.37 

38.10 

«xx 

-74.21 

-59.37 

-44.53 

-29.68 

-14.84 

0 

14.84 

29.68 

44.53 

59.37 

74.21 

The  computed  solutions  are  in  good  agreement  with  the  exact  solution.  For  the 
following  calculations,  we  chose  the  mesh  sizes  k  =  0.001,  h  =  0.05  or  less. 

We  chose  parameter  values  roughly  in  the  range  of  Matta's  experiments, 

reported  in  [8].  (All  the  following  numbers  are  in  CGS  -  units).  We  picked 

L  “  1,  6  =  and  ufl(x)  =  x,  i.e.  we  start  with  a  uniform  filament  of 

M 

length  2  cm  and  radius  1  mm.  We  took  n  =  10,  p  =  1,  — —  =10  and  g  = 

20  (in  Matta's  experiment,  the  drop  falls  into  another  fluid  rather  than  into 
air).  For  the  surface  tension  parameter,  we  chose  the  values  o  =  0,  o  =  1 
and  a  =  10  (the  last  one  is  realistic),  and  for  the  memory  function  we  used 
a  =  0,  a  =  e_t  and  a  =  10e~t.  It  should  be  noted  that  a  =  e~t  increases 
the  shear  viscosity  by  only  10%,  while  a  =  10e-t  doubles  it. 

■ftie  computed  results  for  the  Newtonian  fluid  show  a  pronounced  boundary 
layer  at  x  =  -1.  This  is  illustrated  by  Figure  1,  which  was  obtained  for 
a  •  0,  h  =  0.05  and  k  =  0.001.  As  would  be  expected,  the  effect  of  the  dis¬ 
cretization  is  to  diminish  the  rapid  growth  of  ux  and  the  real  values  for 
ux(-1)  are  even  larger.  This  is  illustrated  by  the  following  table,  which 


compares  values  of  u^-l)  and  uxx(-1)  for  two  different  values  of  h: 


-5- 


u„„(-1) 


Ux(-1) 


Uxx*'1* 


0.7 

7.84 

-25.19 

-0.59 

8.01 

-29.85 

-0.66 

0.8 

12.68 

-66.07 

-0.51 

13.38 

-94.63 

-0.66 

0.9 

20.81 

-162.94 

-0.42 

23.65 

-326.42 

-0.65 

1.0 

33.85 

-362.52 

-0.32 

44.23 

-1201.38 

-0.61 

The  second  column  shows  a  more  rapid  growth  of  and  in  particular  of  uxx. 

u 

XX 

The  quantity  ~r~  should  theoretically  be  independent  of  t,  as  one  can 

ut 

x 

easily  deduce  from  the  equation.  Away  from  x  -  -1,  the  computed  values  are 
almost  identical.  We  believe  it  is  likely  that,  in  fact,  the  values  of 
UxX("D»  u^-1)  tend  to  infinity  in  finite  time,  but  it  would  take  rather 
strong  mesh  refinements  to  verify  this  numerically. 

As  expected,  the  presence  of  surface  tension  accelerates  the  growth  in 
the  boundary  layer.  The  following  table  shows  some  computed  values  of  ux(-1 ) 
and  Uj^-O.S)  (h  -  0.005,  k  «  0.001): 


c  -  0 


a  =  10 


u„(-1)  u_(-0.8) 


44.23 

906.72 

3980.76 


vl,(-1)  u„(-0.8) 


9.62 

16.11 

20.50 


65.63 

2372.01 

8108.24 


77 
9.20 
11.96 
9.43 


u„ ( — 1 )  u„(-0.8) 


1731.19  1.37 
4791.56  0.35 
9085.79  0.16 


The  table  shows  that  the  effect  of  surface  tension  at  c  ■  1  is  relatively 
slight,  but  is  large  at  a  «  10  (a  realistic  value). 


Elasticity  has  the  effect  of  smoothing  out  the  boundary  layer  and 
delaying  the  rapid  growth  of  ux(-1).  The  following  talbe  shows  computed 
values  for  u^-l),  u^-0.8)  at  t  *  1  (h  *  0>005,  k  *  0*001): 


a 

ux(-D 

=  0 

Ux(-0.8) 

a 

ux(-n 

-  1 

ux(-0.8) 

a  = 

ux(-D 

10 

ux(-0.8) 

a(t) 

-  0 

44.23 

9.62 

65.63 

9.20 

1731.19 

1.37 

a(t) 

=>  e-t 

23.38 

10.29 

26.53 

10.21 

68.78 

3.51 

a(t) 

*  lOe-t 

9.63 

8.30 

9.79 

8.36 

13.07 

9.9 

We  can  clearly  see  that  even  small  amounts  of  elasticity  have  a  substanti  1 
stabilizing  influence. 

An  interesting  feature  arises  in  the  case  a(t)  «  e-t,  o  =  1  (Figure 
2).  We  see  a  step  in  the  surface  at  approximately  x  *  -0.1.  In  fact,  the 
computed  values  at  t  *  10  are  ^(-0.2)  =  478.43,  ux(0.0)  «  1.74.  Roughly 
half  the  filament  (x  <  -0,1)  is  stretching,  while  the  other  half  is  actually 
shrinking  and  ultimately  absorbed  into  the  drop.  This  sort  of  feature  does 
not  occur  if  surface  tension  is  absent. 

A  key  to  a  qualitative  understanding  of  the  mechanism  governing  these 
results  may  be  in  results  found  by  Fetrie  [10]  and  Renardy  [5],  [11].  They 
discuss  stretching  of  a  filament  of  uniform  thickness  with  no  inertia.  Hie 
filament  is  stretched  by  a  force  pulling  the  ends.  In  the  Newtonian  case  the 
length  of  the  filament  can  reach  infinity  in  finite  time,  even  for  a  finite 
force.  Hiis  is  because  the  filament  gets  weaker  as  it  gets  thinner,  and  there 
can  be  a  catastrophic  growth.  Elasticity  of  the  fluid  opposes  this,  and  for 
the  model  studied  here,  a  blow-up  is  impossible. 
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In  the  present  situation,  the  stretching  of  the  filament  is  very  inhomo¬ 


geneous  and  this  is  mostly  due  to  the  influence  of  gravity.  The  filament 
"wants"  to  fall,  but  is  prevented  from  doing  so  at  x  =  -1.  Thus,  it 
stretches  most  rapidly  at  this  particular  point.  Our  results  show  a  qualita¬ 
tive  similarity  to  those  mentioned  above:  In  the  Newtonian  case  we  find 
catastrophic  growth  at  the  weakest  point  of  the  filament.  Elasticity 
suppresses  or  at  least  delays  this  effect. 

In  the  previous  calculations,  breaking  of  the  filament  always  took  place 
at  x  =  -1.  In  experiments,  this  is  usually  not  the  case  for  two  reasons: 

1.  The  filament  is  not  really  "attached"  at  the  upper  end.  Instead, 
fluid  is  still  supplied  from  the  nozzle. 

2.  The  assumption  of  a  homogeneous  filament  as  initial  state  is  not 
realistic. 

We  did  some  calculations  with  the  initial  value  uQ(x)  =  1.25x  -  0.25x3, 
corresponding  to  a  filament  that  is  thinner  in  the  middle.  Otherwise  we  used 
the  same  parameters  as  above.  The  following  table  shows  the  values  of  u 
at  t  »  2  for  a  Newtonian  filament  and  a  =  0  and  10  (the  step  sizes  are 
again  h  =  0.005,  k  =  0.001): 


X 

u(x,t=2) ,o=0 

u(x,t=2) ,0=10 

-1 

-1.00 

-1.00 

-0.8 

-0.21 

-0.98 

-0.6 

7.22 

-0.95 

-0.4 

17.97 

-0.92 

-0.2 

26.30 

-0.85 

0 

31.75 

29.10 

0.2 

34.86 

29.14 

0.4 

36.50 

29.17 

0.6 

37.37 

29.19 

0.8 

37.83 

29.21 

1.0 

38.06 

29.23 

The  second  column  shows  that  the  filament  breaks  somewhere  between  x  **  -0.2 
and  x  *  0. 
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